Magnetic field induced localization in 2D topological insulators 
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Localization of the helical edge states in quantum spin Hall insulators requires breaking time 
reversal invariance. In experiments this is naturally implemented by applying a weak magnetic field 
B. We propse a model based on scattering theory that describes the localization of helical edge states 
due to coupling to random magnetic fluxes. We find that the localization length is proportional to 
B~ 2 when B is small, and saturates to a constant when B is sufficiently large. We estimate especially 
the localization length for the HgTe/CdTe quantum wells with known experimental parameters. 
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The prediction and discovery of quantum spin Hall insula- 
tprs (QSHIs) [TJ[TJ[51|5] nas opened a door to an unexpected 
Category of topological phases in condensed matter [SHE], 
C^Lnd revealed a new route to investigations of edge /boundary- 
T^tate physics [9HTT]. Although the prototypes of QSHIs [HQ] 
^^,re mainly based on two copies of quantum Hall insulators, 
—which have been investigated for more than three decades 
, " 1,1 v| [T5] , it was soon realized that the fundamental impor- 
tance of time reversal invariance (TRI) distinguishes the two 
1 systems in a profound way [2j. Indeed, QSHIs, unlike the Z- 
'"^lassified quantum Hall insulators [2], belong to a class of 
r^ wo-dimensional time-reversal-invariant Z2 topological insu- 
lators [2J. The defining feature of QSHIs, as its name sug- 
gests, is a pair of helical edge states that persist in the bulk 
gisulating gap of the system [TJ [TJ EJ [9] . 
| » The topological power of QSHIs lies precisely in the ro- 
bustness of the helical edge states against generic pertur- 
bations due to unavoidable disorder in every experimental 
isetup, unless TRI is broken. In the presence of both TRI 
breaking and disorder, the helical edge states will be lo- 
calized, and the general framework of Anderson's localiza- 
tion theory applies [E]. Nevertheless, the localization of 
Lhc helical edge states distinguishes itself from conventional 
o ne-dimensional localization when the focus is placed on the 
^rucial role TRI plays in the problem. This point becomes 
("^especially relevant as TRI can be broken continuously, for 
^instance, by turning on a magnetic field gradually. Indeed, 
^he sensibility of transport though helical edge states to weak 
magnetic field has been demonstrated experimentally in the 
^measurement of magneto-conductance in topologically non- 
^jrivial HgTe / CdTe quantum wells [SJ [3] . Related theoretical 
analyses have been carried out that include the interplay be- 
• bween TRI-breaking and disorder, but mainly consider mag- 
• £ijetic impurities [17[ 118] , or bulk random potential combined 
S^rith magnetic field [TH]- A transparent edge theory that 
SJbcuses on the magnetic-field-dependent localization of the 
^nelical edge states, however, is still missing. 

In this paper we propose a model that explicitly addresses 
the question on how the localization of helical edge states 
occurs as a weak magnetic field is gradually turned on. Our 
model is based on the scattering theory of edge states in the 
presence of generic edge disorder. In particular we consider 
the existence of alternative paths for the edge states due 
to, e.g., constrictions formed at the rough edges of realistic 
samples (see Fig. [TJ, which further allows for loops of the 
helical edge states. The magnetic field penetrating though 
these loops results in broken TRI that is experienced by the 
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FIG. 1: (a) Helical edge states in a disordered QSHI in a uni- 
form magnetic field. Occasional occurrences of constrictions along 
the edge lead to Fabry-Perot-type loops where Aharonov-Bohm 
phases due to magnetic fluxes can accumulate, (b) The scattering 
of the helical edges by one of these loops, described by a scattering 
matrix S, can be divided into two parts: the scattering between 
two pairs of helical edge states (S), and the propagation of one 
of these pairs around the loop (S$). (c) Three types of scatter- 
ing probabilities, T, R and T s , that are relevant to the scattering 
between two pairs of helical edge states. 



helical edge states in the form of finite random fluxes. We 
show that these random fluxes necessarily lead to localiza- 
tion of the helical edge states with universal behaviors in two 
regimes: immediately after the magnetic field is turned on, 
the localization length becomes finite and decreases as B~ 2 : 
when the standard deviation of the random fluxes (propor- 
tional to B) is comparable to or larger than one magnetic 
flux quantum, the localization length saturates to a constant. 
In-between these two regimes, damped oscillations of the lo- 
calization length may arise depending on the distribution of 
the random fluxes. Our results provide a clear illustration of 
the symmetry-breaking-induced localization in QSHIs. 

We start to introduce our model by considering one of its 
possible realizations, depicted in Fig. [T] The edge roughness 
of a realistic QSHI sample may lead to occasional constric- 
tions (at one edge) where the helical edge states can either 
tunnel across or pass around. As a consequence, loops can 
form and attach to the propagating path of the helical edge 
states. When a perpendicular magnetic field is applied to 
the sample [25] , each of these loops acts as a magnetic flux 
impurity, to be distinguished from usual magnetic impuri- 
ties. Individually such an impurity works like a Fabry-Perot 
scatterer (see Fig. [IJd), where the scattering probability am- 
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plitudes depend on the Aharonov-Bohm (AB) phase, owing 
to the magnetic flux, acquired by electrons circling around 
the loop. The collective action of a random distribution of 
these scatterers causes localization of the helical edge states 
with an explicit dependence on the magnetic field. The main 
part of this problem can be tackled consistently by scattering 
theory, as we now show. 

To analyze the scattering of the helical edge states by a 
single magnetic flux impurity (see Fig. [T|, we divide the full 
scattering process into two effective parts: the local scatter- 
ing between two pairs of helical states, and the free propa- 
gation of one pair of helical states that closes the loops. For 
simplicity, we assume that the first part does not depend on 
magnetic field, hence respects local TRI, while the second 
part contains the entire information about the magnetic flux 
by means of AB phases that enter the propagators. 

Owing to the local TRI, the scattering between two pairs 
of helical states, described by a 4 x 4 scattering matrix S, 
has the following constraint: 
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where is the time-reversal operator. We choose a specific 
basis ordered as (Ri, L\, L%, R2), where Ri (Li) stands for 
the right (left) mover of the i-th Kramers pair (i — 1,2), 
such that the time-reversal operator reads 
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with o v the Pauli matrix and k the complex-conjugate op- 
erator. Consequently, the scattering matrix S satisfying Eq. 
necessarily has the following form: 



S = 



fh r' s'\ 

ti -s r 

r -s' t 2 
V s r' t 2 ) 
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Here U stands for the direct transmission for Kramers pair 
i; r (r 1 ) stands for the reflection from a right (left) mover 
to a left (right) mover; s and s' represent the transmission 
by switching to another Kramers pair. Importantly, zeros 
in S signify the absence of direct back-scattering within one 
Kramers pair due to TRI. Taking into account the unitar- 
ity of the scattering matrix, the parametrization of S can 
be further simplified as (up to an unimportant global phase 
factor): h = -t* 2 =t,r' = r* , s' = s* , and T + R + T s = 1, 
where T = \t\ 2 , R = \r\ 2 and T s = \s\ 2 . 

The free propagation of Kramers pair 2 leads to a 2 x 2 
scattering matrix S$ that is diagonal in the basis (£2,-^2), 
given by 
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with ip the dynamical phase and <f> the AB phase (equal to 
the total flux enclosed by the loop in units of h/2ne). 

Combining the two parts above, we find the final scattering 
matrix S for Kramers pair 1, in the basis (Ri, L\), to be 



RZ+ + T S Z_ 
-rsAZ 



-(rs)*AZ 
\t\+RZ_+T s Z, 
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where 



Z+ = 



a i(ip±<p) 



1 + |f|e*(v d 



AZ = Z, 



Z_ 



\t\ cos 4> + cosh(ln \t\ + iip) 



(6) 
(7) 



and the phase of t* has been absorbed into the dynamical 
phase if. The back-scattering probability can be obtained 
immediately: 
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Evidently, for the helical edge states to be back-scattered 
with finite probability, two conditions must be fulfilled. 
First, it is necessary that both R and T s are finite. If one of 
these two tunneling probabilities is zero, the system essen- 
tially reduces to two decoupled copies of quantum Hall edge 
states, and back-scattering is known to be prohibited for ei- 
ther copy [20 . Only when both tunneling processes (i? and 
T s ) are allowed, the system belongs truly to the Z2-classified 
symmetry class where TRI plays a central role. It follows 
that the second necessary condition for back-scattering to 
occur is to break the global TRI by having 0^0 mod n. 
The cooperation of these two conditions clearly illustrates 
the underlying protection mechanism, from the scattering 
point of view, for the helical edge states of QSHIs. 

At disordered sample edges, the magnetic flux impurities 
will occur randomly, and the helical edge states will even- 
tually be localized as a consequence of finite back-scattering 
probabilities for individual scatterers. Here, we assume not 
only that the variables (including phases and scattering am- 
plitudes) for each individual scatterer are random, but also 
that different scatterers are completely independent such 
that the relative scattering phases for two consecutive scat- 
terers are uniformly distributed. The localization length of 
the helical edge states in this scenario can be extracted from 
the appropriate scaling variable lnT, where T = 1 — TZ 
is the total transmission probability for the effective one- 
dimensional system |21| . 

The total transmission probability Tn through N scatter- 
ers can be calculated by multiplying the transfer matrices 
that relate the scattering amplitudes on the right side of 
each individual scatterer (labeled i, i = 1, N) to those on 
the left. A general transfer matrix reads 



M, 



1 
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Aj Pi 

Pi A i 
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where Tj corresponds to the transmission amplitudes (diag- 
onal entries in S) for the «-th scatterer, pi corresponds to 
the reflection amplitudes (off-diagonal entries in S), and Aj 
is a phase factor that is independent for each scatterer. Note 
that the dynamical phase for the free propagation of the he- 
lical states in-between two consecutive scatterers (i and i + 
say) can be obviously incorporated into the above transfer 
matrix while preserving its general form. Note also that 
multiplications of the transfer matrices preserve the general 
form as well. We will put an overhead tilde to distinguish 
the amplitudes involving i consecutive scatterers from those 
involving only a single (i-th) scatterer. Then T/v is given by 
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with e lB = \n-i^nPn-iP*n I\Pn-iPn\- Our assumption 
of independent scatterers implies that 9 is uniformly dis- 
tributed in [0, 2ir). It follows that the mean of the scaling 
variable (InT) becomes simply additive |21| 



(In Tat) = (ln71v-i) + (lnM 2 ), 



(11) 



where (ln|l + \pN-iPN\ el9 \ ) has vanished after averaging 
over 9. 

The inverse localization length 7 = \/i\ oc is denned in 
terms of the scaling variable as: 
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-n(ln(l-ft)), 



(12) 



where n is the linear density of the scatterers. The final av- 
erage (ln(l — 1Z)} is over certain distributions of independent 
variables including <p, (f> and scattering amplitudes that ap- 
pear in Eq. Q for a single scatterer. We are interested in 
the weak back-scattering case for each individual scatterer 
(IK. 1), thus 



7 ~ n(TZ). 



(13) 



The average in terms of the arbitrary dynamical phase ip can 
be carried out exactly, and yields 



RT S 1+T 



R + T s T 



sin 



(1-T) 2 
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By further using the fact that <fr = BA(2ire/h), where the 
magnetic field B is taken to be uniform and A is the area 
enclosed by the helical loops, we will only need to average 
over distributions of the scattering amplitudes and the area 
A in order to estimate the localization length l\ oc = 1/7. 



One immediate consequence of Eq. ( 14 1 is that the localiza- 



tion length is magnetic field symmetric, which is certainly 
expected. 

An important regime that we are particularly interested in 
is the weak magnetic field regime, where BA -C h/e with A 
the mean of A. In this regime Eq. ( 14 1 becomes (assuming 
T is not too close to 1) 



7 = aB 2 

„ /27re\ 2 / RT S 

with a = 4n[ ( 

V h ) \(R + T S 



(15) 

(1 + T))(A 2 ). (16) 



Manifestly a is a constant factor for given distributions of 
scattering amplitudes and A. The independence of the 
inverse localization length here is a universal result of our 
model in the sense that it does not depend on the specific 
distributions of variables for individual scatterers. It im- 
plies that the localization length of the helical edge states 
is finite at weak magnetic field and diverges only as 1/B 2 
when the magnetic field is vanishing. Furthermore, the 
low-temperature magneto-conductance of a QSHI should 
also vary as B 2 in the weak magnetic field limit, that is, 
5G(B) = G(0) - G(B) oc B 2 . This contrasts our result with 
the linear magneto-conductance behavior previously found 
on a lattice model with bulk impurity potentials |19| . 

Another interesting regime is when the magnetic field is 
strong enough such that both BA and Bo a (with a a = 
\J ((A — A) 2 } being the standard deviation) are comparable 




FIG. 2: T(B) for different values of T, where A is assumed to 
take a normal distribution with a a/ 'A = 0.3. 



or larger than a flux quantum h/e. In this regime we can 
approximate the average in terms of A as an average over a 
uniform distribution of </>, which yields 



7 sa t = 2n 



RT S 
R + T s 



(17) 



This especially simple result again shows a universal behav- 
ior in our model-the inverse localization length saturates 
at relatively strong magnetic field irrespective of the spe- 
cific distributions of scattering variables. However the actual 
value of 7 sat certainly depends on the distributions of R and 
T s . Moreover, the above formula re-emphasizes the impor- 
tance of allowing both tunneling processes represented by R 
and T s to evoke a true protection mechanism due to TRI. 

In-between the two regimes discussed above, we need 
to consider the specific distributions of the variables in 
Eq. (14). Let us first focus on the behavior of 7 by 
assuming that A has a Gaussian distribution character- 
ized by the mean A and the standard deviation a a- It 
is instructive in this case to look at the function T(B) = 
[(1 + T)/2T] (sin 2 0/[sin 2 (j>+ (1 - T) 2 /4T]) A with the aver- 
age only taken in terms of A. T(B) has been defined such 
that it saturates to the value 1 at sufficiently strong magnetic 
field. In Fig. [2] we plot the numerically obtained T(B) for 
various T and fixed A and a a ■ Right after the magnetic field 
is turned on, T(B) shows a quadratic increase irrespective of 
the assumed T or the distribution of A. Before T(B) satu- 
rates, it undergoes damped oscillations when a a/ A is small. 
These oscillations are due to the collective AB effect for the 
helical loops in our model: the loops enclosing similar area 
lead to AB oscillations of similar period; they contribute co- 
herently to the back-scattering of the helical edge states; the 
magnetic field dependence of the total transmission is thus 
shaped by the AB effect at individual scatterers when Bua 
is significantly smaller than <fio = h/e. The period of the os- 
cillations is roughly 0o/2A, where the factor 1/2 is obviously 
a consequence of T (and hence 7) only depending on sin 2 (f>. 
[27 The overall amplitude of the oscillations is suppressed 
for large T and enhanced for small T. The reason is intu- 
itively clear: the more the helical edge states are scattered 
into the loops, the more pronounced the resulting AB effect. 

Now we address the issue of the scattering amplitudes 
which have only been assumed to be phenomenological pa- 
rameters in the scattering matrix S so far. To this end we 
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FIG. 3: A schematic view of a constriction where two pairs of 
helical edge states (indicated by the linear bands) are coupled 
(indicated by the mixed bands) in a region of length L and sepa- 
ration W. 



investigate a constriction depicted in Fig. [3] which is de- 
scribed by the following effective Hamiltonian: 
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hv F k x J 



(18) 



where vf is the Fermi velocity for the helical edge states, 
m(x) and 5 (x) represent x-dependent coupling between the 
edge states, and the basis is ordered as (Ri, L\, L 2 , Rz)- The 
above Hamiltonian manifestly respects TRI: H — OHO -1 . 
This Hamiltonian can be derived from microscopic mod- 
els such as the Bernevig-Hughes-Zhang (BHZ) model for 
HgTe/CgTe quantum wells [THU Hfl |55]. For simplicity we 
take m(x) = mwd{x)8(L — x) and S(x) = SwO(x)9(L — x) 
with 9(x) the Heaviside step function and mw and Sw two 
constants determined by the constriction separation W. In 
the case of HgTe/CgTe quantum wells, a nonvanishing S term 
owes its existence to the presence of bulk-inversion asymme- 
try m ig. 

The scattering amplitudes for this constriction, corre- 
sponding to S in Eq. (j3j, can be easily derived (see sup- 
plementary materials): 



t = icos(5wL/hvF)/Ci 
s = s'm(S w L/hv F )/(, 
r = mw s'm(qL) / hvpq(, 



(19) 
(20) 
(21) 



where q = ■sj E 2 — m 2 ^ /Hvp can be either real or imaginary 
depending on the energy E. and £ = \(E/HvfQ) sin(qL) + 
icos(qL)\ is a normalization factor. Clearly r vanishes when 
mw — 0, which shows the fact that m couples R\ (L±) to 
L2 (R2)', s vanishes when Sw — 0, which shows the fact 
that S couples Ri (Li) to R 2 (L 2 ). For low energy (\E\ < 
mw) scattering states, R/(T + T s ) ~ sinh 2 (mw L / Kv f) if 
\E\ < m w , and R/(T + T S ) ~ (mwL/hv F ) 2 if \E\ ~ m w , 
meaning that the reflection dominates as long as L is large 
compared with hvp/mw- In this regime, Eq. (12 1 reduces 
to 7 ~ 4n(T s )(sin 2 4>). 

More generally the average with respect to the scattering 
amplitudes has to be performed numerically by taking cer- 
tain distributions of W and L (at a certain energy E). The 
advantage of this change of variables is that W and L, un- 
like the scattering amplitudes, are in principle independent 
to each other. In total, this leads to three independent ge- 
ometric variables, W, L and A, that remain to be averaged 




FIG. 4: Inverse localization length 7 as a function of magnetic 
field B, for different distributions of loop area A. It shows a B 2 - 
increase in the weak B limit (inset) and a saturation at relatively 
high field. In-between, damped oscillations may occur if a a/ A < 
1. The saturation value 7 sat is sample-dependent but has an order 
of magnitude O.Oln (with n the linear density of the scatterers) 
in our realistic estimations. 



on in our final evaluation of 7 as a function of magnetic field 
B (we will not make any assumption on the density n of 
scatterers and leave it as a parameter). After carrying out 
these averages numerically (see supplementary materials for 
details), we show the typical results in Fig. [4] The qual- 
itative behavior of j(B) in Fig. [4] is essentially the same 
as that of T(B) in Fig. [2] except that 7(B) is shown for 
various ratios a a/ A whereas the scattering amplitudes have 
been averaged out. The universal features which we can ob- 
serve are that 7 increases as B 2 at weak magnetic field and 
saturates at sufficiently strong magnetic field. Despite the 
fact that the exact value of 7 sat depends on the energy E 
and the distributions of W and L, the order of magnitude 
°f 7sat turns out to be consistently O.Oln for all cases with 
realistic considerations (see supplementary materials). We 
also observe in the intermediate regime damped oscillations 
of 7 which are pronounced if a a/ A is small but suppressed 
as long as a a/ A is close to or larger than 1. We point out 
here that 7(B) has a local minimum/maximum, hence the 
localization length has a local maximum/minimum, when- 
ever B is roughly an integer/half-integer multiple of <j)Q/2A 
- this is where the TRI is maximally preserved/broken. 

To summarize, we have investigated a simple yet illumi- 
nating model that demonstrates a magnetic-field-dependent 
localization of the helical edge states in quantum spin Hall 
insulators. We have identified universal, sample-independent 
features, as well as an interesting but sample-specific feature 
in this model. With known parameters for the HgTe/CgTe 
quantum wells, we have also estimated quantitatively the lo- 
calization length. Both the qualitative and the quantitative 
results can be examined by experiments. 
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I 

SUPPLEMENTARY MATERIAL 



DERIVATION OF EQS. (19-21) 

In this section we derive the scattering amplitudes, given by Eqs. (19-21) in the main text, for the constriction described 
by the Hamiltonian (18) in the main text. By assuming m(x) = mw9(x)9(L — x) and S(x) = S\y9(x)9(L — x) with 9(x) the 
Heaviside step function and mw and Sw two constants, we divide the constriction into three regions: x<0,0<x<L and 
x > L. The scattering amplitudes are obtained by matching energy eigenstate wavefunctions for adjacent regions. 

The energy eigenstates in the regions x < and x > L are trivial: 
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C3e -ik x 

\ c 4 e lk ° x 



(22) 



where ko = E/hvp. In the region < x < L, the energy eigenstate is given by 
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(23) 



v -i / 



where = a,rccos(E/m w ), k\ = —k 2 = Sw/hvp + q and fc 3 = — fc 4 = ~8mjhv F + q with q = \fE 2 — m^/hvp = 
i(mw /Twf) sin 8. Note that both 9 and q can be complex depending on the energy E. Note also that the basis for the above 
wavefunctions is ordered as (R\, L\, L 2 , R2) (see Fig. 3 of the main text). 

To start with we assume that the only incoming state is from channel R\, that is, a\ = 1 and a 4 = c 2 = c 3 = 0. Then we 
need to match the wavefunctions such that 



* b (.T = 0) = * a (a; = 0) = 



^ b (x = L) = ^ c {x = L) = 
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FIG. 5: Schematic setup for the simulations. 
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a (nm) 


364.5 


-686 


-7.5 


-512 


-10 
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TABLE I: Parameters for the simulations. Note that the parameter C, which corresponds to an overall constant energy shift, takes 
the value so that the edge bands cross at E — 0. 

The above equations fix the values of bi (i = 1,2,3,4), and hence the values of 02, 03, c\ and C4. In particular, ai can be 
identified as the backscattering amplitude and we find a 2 = identically; a 3 can be identified as r; ci can be identified as t; 
C4 can be identified as s. We find: 

t = i cos(SwL/Hvf) sin 0/ sin(qL — 9), (26) 
s = s\n(5w L / hv f) sin6> / ' sin{qL — 0), (27) 
r = -ism(qL)/shx(qL-9), (28) 

which are precisely Eqs. (19-21) in the main text after removing an unimportant common phase factor. By assuming 
differently the incoming states, we can construct the full scattering matrix for the constriction. Both the unitarity of the 
scattering matrix and the symmetry relations between the matrix elements can be easily checked. 



NUMERICAL SIMULATIONS TO EXTRACT THE DISTRIBUTIONS OF SCATTERING PROBABILITIES 

In this section we extract from numerical simulations the distributions of the scattering probabilities for the constriction 
illustrated in Fig. 3 of the main text. In our simulations we employ a six-terminal setup as shown in Fig. [5] This setup is 
equivalent to a Hall-bar setup. We define a point contact with cosine profiles to simulate the effect of the constriction. The 
point contact has two parameters: its length L and its separation W. The depleted regions are defined by a sufficiently high 
on-site potential (compared with the band width). 

The model Hamiltonian for the central region of the setup is the tight-binding Hamiltonian corresponding to the Bernevig- 
Hughes-Zhang (BHZ) model Q: 

with h(k) = (C- Dk 2 ) + Ak x a x - AkyO y + (M - Bk 2 )a z (30) 

and (Tj (i — x, y, z) are Pauli matrices. The block-off-diagonal term, proportional to A, is a spin-orbit interaction term due to 
the bulk inversion asymmetry in HgTe/CdTe quantum wells |5J[3]. We list in Table [i] the experimentally obtained parameters 
for the above Hamiltonian, as well as the lattice spacing a adopted to discretize this Hamiltonian. We will measure length 
in units of a hereafter. With these parameters we estimate the Fermi velocity for the helical edge states to be vf — 3.8 x 10 5 
m/s. 

The scattering probabilities for the constriction that are needed in the main paper can be identified in the current setup 
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FIG. 6: Typical dependences of scattering probabilities on W and L. For the left panel, L = 30a; for the right panel, W = 10a. 




FIG. 7: Histograms showing the distributions of the scattering probabilities. 



as follows: 

T = T 21 ,T S =T 3U R = T 41 , (31) 

where Tji is the transmission probability from contact i to contact j. These transmission probabilities can be calculated by 
using the standard Green's function technique [1]. We also check the sum rule T sum = + T31 + T41 = 1 to ensure the 
validity of our results. 

For the distributions of the scattering probabilities we need to assume reasonable distributions for L and W. Without 
sufficient knowledge from experiments we make the following assumptions: L is uniformly distributed in the range (0, 21 go) 
with the spin-orbit length lso — hvp/A ~ 30a; W is uniformly distributed in the range (0, 2£) with £ the penetration depth 
of the edge states. £ is energy dependent and its order of magnitude is given by Hvf/M ~ 5a. We will fix the energy at 
E = 0.5 meV for the results presented below, where £ ~ 10a. We have checked the robustness of our results with other values 
of energy and/or other reasonable distributions (e.g. Gaussian distribution) of L and W. 

As results of our simulation, we show the typical dependences of the scattering probabilities on L and W in Fig. [6] and 
the histograms for the distributions of the scattering probabilities in Fig. [7J To obtain our final result presented in Fig. 4 
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of the main text, we simply generate samples with randomly chosen L and W according to our assumptions; no accurate 
distribution functions for the scattering probabilities are actually needed for our purpose. 
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